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ABSTRACT 



We study the three dimensional fundamental-adjoint SU(2) lattice gauge 
theory at finite temperature by Monte Carlo simulations. We find that the 
finite temperature deconfinement phase transition line joins the first order 
bulk phase transition line at its endpoint. Moreover, across the bulk tran- 
sition line, the Polyakov loop undergoes a discontinuous jump implying the 
existence of both confining and deconfining phases on its two sides. Implica- 
tions for universality and the nature of the confining-deconfining transition 
are discussed. 
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1. INTRODUCTION 



Quantum field theories are defined with an ultraviolet cut-off. The choice 
of this regulator is highly non-unique. In perturbation theory the renormal- 
izability criterion ensures that the physical answers are independent of the 
cut-off parameter used. In fact, different choices of the regulators only change 
the scale in the theory, leaving physics unchanged. 

A non-perturbative way to regulate a theory is to define it on a lattice, 
which has a built-in ultraviolet cut-off (the lattice spacing a). A regulated 
theory on the lattice corresponding to a particular theory in the continuum is 
non-unique due to the large freedom of choice of lattice action and the lattice 
type. This non-uniqueness is expected to be a lattice artifact and should dis- 
appear as the lattice spacing a tends to zero (or the lattice correlation length 
— > oo), resuling in a universal continuum limit. The Wilson formulation [|TJ 
of lattice gauge theory is a particular choice of lattice action from a much 
more general class, given by 

S = Y.Y.I 3 R ReTr R U P • (!) 

P R 

Here Up, the plaquette action, is an ordered product of the four group el- 
ements Uj (e SU(N)) corresponding to the four links of an elementary 
plaquette. The sum over R extends over all the representations of the gauge 
group and the sum over P runs over all plaquettes. Ttr denotes the trace 
over colour space in the representation R. 

A variety of different forms of lattice actions have been employed in the 
past to check the scaling behaviour and the universality of the SU (2) gauge 
theories. Dimensionless ratios of physical quantities, such as Mc/y/o 1 ^, or 
T c /-y/a|J, where Mq is the lightest glueball mass, T c is the deconfinement 
temperature and a is the string tension, have been computed for the SU(2) 
lattice gauge theory and shown to be almost independent of the lattice ac- 
tions used. Even more stringently, functions of such dimensionless ratios, 
e.g., the energy density as a function of T/T c , have been successfully used to 
demonstrate [§] universality. While these results have boosted our confidence 
in the existence of a non-perturbative continuum limit of the SU(2) gauge 
theory, it is still not clear whether more qualitative aspects of the theory, 
such as the order of the deconfinement phase transition, also exhibit this 
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universality. We investigate the confining- deconfining behaviour in the ex- 
tended coupling plane of a particular class of theories belonging to eq. (jl|) 
and compare it with that of the usual Wilson action. The organization of 
this paper is as follows: In sect. 2 we introduce the model and the motiva- 
tion to study it at non-zero temperatures. In sect. 3 we give a very brief 
review of the observables computed and their critical behaviour. In sect. 4 
we discuss the data and its analysis. Sect. 5 consists of the observations and 
their implications. 



2. THE MODEL 



A well studied model belonging to the above class of SU (2) gauge theory 
is described by the action 



p 



i 



Tr F U P )+p A l--Tr A U P 



(2) 



Here F and A denote the fundamental and adjoint representations respec- 
tively. Comparing the naive classical continuum limit of eq. (0) with the 
standard SU(2) Yang-Mills action, one obtains 



J_ = /3 Pa 
9l 4 + 3 



(3) 



Here g u is the bare coupling constant of the continuum theory. One can also 
characterize the model with a pair of couplings (g u ,6) with tan# = Pa/P- 
In the non-perturbative continuum limit (i.e. a — > 0) each of these theories, 
characterized by a given 9, flow to the same critical fixed point, = 0. The 
corresponding asymptotic scaling relation || is 



A(9) 



exp 



fort 



-0i 



(4) 



where 



log 



A(0) 5vr 2 6tan0 



A(0) 11 (3 + 8tan( 



(5) 
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Here (3q and (3\ are the usual first two coefficients of the (3 function for the 
SU (2) gauge theory. 

This model has been extensively studied after Bhanot and Creutz || 
found it to have a rich phase structure (fig.l). At (3a — 0, it reduces to the 
standard Wilson action. Along the (3 = axis it describes the SO (3) model 
which has a first order phase transition at (3 C A ~ 2.5. At (3a = oo it describes 
the Z 2 lattice gauge theory again with a first order phase transition at (3 C = 
\tn{\ + y/2) ~ 0.44 0. Bhanot and Creutz found that these first order 
transitions extend into the ((3, (3 a) plane, ending at an apparent critical point 
located at (1.5,0.9). These transition lines are shown in fig.l by continuous 
lines. The same phase diagram was also later produced qualitatively by 
mean field theory techniques ||. The A-ratio was found to be smaller than 
its asymptotic value given by eq. (0) by a factor of ~3.5 in a Monte Carlo 
determination from the ratio of string tension evaluated on (3a =constant 
lines around (1.5, 0.9) and (2.2, 0.0) in the ((3, /?4)-plane. It was later 
shown 0| that the 9 dependance of this factor further differs rather strongly 
from eq. ([5]) on the (3a = 0.9 line. Thus in the vicinity of the endpoint of 
the first order line, the asymptotic scaling relation is known to be violated 
strongly. Ref. M also found that the string tension does not go through zero 
in this region, suggesting the endpoint to be a higher order phase transition, 
although the location of endpoint was obtained in Ref. || by extrapolating 
the discontinuity in average action to zero. 

The motivation for Bhanot and Creutz to study this model by Monte 
Carlo simulation was to show that non-abelian lattice gauge theories need 
not be free of singularities to maintain confinement in the asymptotic scal- 
ing regime. As fig. 1 shows, one can go smoothly around the line of first 
order singularities without affecting confinement in the asymptotic regime. 
However, this is perhaps an incomplete picture of the phase diagram of the 
model. Since the numerical results, which established the phase diagram 
were obtained on small, N*, lattices for N a = 5-7, one must worry about a 
possible deconfinement phase transition on these lattices at sufficiently small 
coupling g u . As all finite lattices are necessarily at finite temperature as well, 
the confining-deconfining transition line must exist in the extended plane on 
all lattices and one needs to understand its behaviour on larger lattices to 
draw conclusions about the zero temperature phase diagram of the theory 
and its possible implications to the continuum limit of the confining phase. 
The deconfinement phase transition has recently been extensively studied at 
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(3a = |T0| , H| . We follow this transition in the (3a direction and find that 



it merges with the bulk phase transition line found by Bhanot and Creutz. 
The phase diagram therefore splits into seperated confining and deconfining 
phases. 



3. THE OBSERVABLES AND THEIR CRITICAL BEHAVIOUR 



On a finite lattice of size x N T the volume and temperature are given 
bypl 



V = (iV» 3 and T = — — . (6) 

1 \ q-CL 

Here a is the lattice spacing governed by the bare coupling g u for sufficiently 
small a, as shown in eqns. (§-[5]). To avoid big finite spatial volume effects, one 
chooses N a > N T for finite temperature simulations. One basic observable 
which we studied is the average plaquette energy (P), where 

l EpTr F U P 

2 QN^N T y } 

Due to periodic boundary conditions in the timelike direction, the SU(2) 
lattice gauge theories described by eq. ([D and in particular eq. @ have a 
Z 2 invar iance at finite temperature corresponding to the center of the group. 
Under this 

U (n,r ) — > zU (n,r ) Vn, r : fixed , and z E Z 2 . (8) 

Here Uo(n, r) is the timelike link at the lattice site (ft, r). This transformation 
leaves the action invariant but the Polyakov loop, defined by, 

I n t 

L(n) = -Trl[U (n,T), (9) 

Z r=l 

changes: 

L -> zL . (10) 
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A nonvanishing value for (L) signals a spontaneous break-down of the global 
Zi symmetry. The thermal expectation value of the Polyakov loop or its 
average value L = Y^L(fi) can also be shown to be the order parameter 

n 

for the deconfmement phase transition since it is a measure of the free energy 
of an isolated free quark ]12, |I~4]j . 

It is been argued that the effective theory for this order parameter, ob- 
tained by integrating out all other degrees of freedom, is in the same uni- 



versality class as the Ising model in 3-dimensions [|1^, [13| . Thus, if the 
deconfmement phase transition is of second order, then it will have the same 
critical exponents as those of the Ising model, provided the universality class 
does not have another second order phase transtion. On an infinite lattice, 
these exponents are /3, 7 and v corresponding to the order parameter itself, 
its corresponding susceptibility and the correlation length: 



(11) 
(12) 

(13) 



For the 3-d Ising model (3 ~ 0.325, 7 ~ 1.24 and v 0.63. For the finite 
temperature SU(2) gauge theory, one obtains these exponents from Monte 
Carlo simulations by simply fitting the order parameter [It], |TJ, [T5| or by 
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using the finite size scaling theory Hl6| for the susceptibilty. According to the 
latter, the susceptibility on a lattice of spatial extent N a is expected to grow 
like 

X^K , (14) 

where u = j/u = 1.97 according to the universality prediction above. If 
the phase transition were to be of first order instead, then one expects the 
exponent uj = 3, corresponding to the dimensionality of the space [jDJ]. In ad- 
dition, of course, the order parameter is expected to exhibit a sharp, or even 
discontinuous, jump and the corresponding probability distribution should 
show a double peak structure. For (3a = 0, the universality prediction 
was verified by Monte Carlo simulation by Engels et. al. |l0fl , who found 
uj = 1.93 ± 0.03. 
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4. DATA AND ANALAYSIS 



Our Monte Carlo simulations were done on 8 3 x 4 and 10 3 x 4 lattices, 
using the Metropolis algorithm. Around (3 a — 11 where the nature of phase 
transition changes from 2nd order to 1st order, we also simulated the model 
on a 12 3 x 4 lattice. The Ferrenberg-Swendsen |H| technique was used to 
extrapolate the data at neighbouring couplings to locate the critical point. 
Once located, much longer runs were performed at these points. (200,000, 
150,000 and 50,000 sweeps on 8 3 x 4, 10 3 x 4 and 12 3 x 4 lattices respectively). 
Table 1 shows the values of the couplings ((3, (3a) where simulations were 
performed along with the critical couplings (/3 C ) determined from the peak 
in susceptibility using Ferrenberg- Swendsen method. The last entry in the 
table (1) is the value of finite size scaling exponent (u) calculated from the 
values of the susceptibility peak on 8 3 x 4 and 10 3 x 4. At ^=0.9 and 1.1 
the second value of u is calculated from the susceptibility peaks on 8 3 x 4 
and 12 3 x 4 lattices. 

As the expectation value of the Polyakov loop is always zero on a finite 
lattice, we measure the absolute value of L after taking its lattice average. 
The values of average Polyakov loop (\L\), extrapolated by the Ferrenberg- 
Swendsen technique in the neighborhood of critical points are plotted in 
fig. 2-a,b,c,d for all (3a along with the actual values determined from the 
individual runs. The error bars shown were determined from binning. As 
seen from fig. 2, the slope of the Polyakov loop curve keeps increasing with 
increasing value of (3a- All the corresponding values of susceptibilities, 

X = Nl{{L*)-(\L\)% (15) 

are plotted in fig. 3-a,b,c,d. We have also shown the first order and second 
order predictions for the susceptibility peaks on the 10 3 x 4 and 12 3 x 4 
lattices based on the measured susceptibility peak on the 8 3 x 4 lattice and 
uj = 3.0 and 1.97 respectively. As is clear from the susceptibility plots, the 
universality hypothesis is well supported by the data for all the (3a including 
and up to (3a = 0.9. Instead of using the critical exponents above, one can 
try to obtain them from the data by simply comparing the peak heights. The 
corresponding exponents from the data for susceptibility are given in table 
1. They exhibit an excellent agreement with the universality prediction of 
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uj = 1.97 for (3a < 0.9, confirming the second order nature of the phase 
transition. 

At Pa — 11 there is evidence that the transition has become first order. 
The histograms of \L\ (fig. 4) on N a = 10, 12 show a definite two peak 
structure with the dips between the two peaks increasing with N a . The 
value of the susceptibility peak, however, lies below the first order predicted 
values but definitely outside the second order predicted region also. This is 
unlike the clear second order nature of the transition up to (3a = 0.9. It is 
possible that at (3a = 11, the corrections to the leading scaling prediction 
of eq. ([14]) are large and thus larger lattices are needed to ascertain the true 
nature of the phase transition. It has been found in simulations of Potts 
models with first order transitions, that much larger lattices are needed to 
see the correct finite size scaling exponents than are needed to see a two peak 
histogram JTJJ. To verify the onset of the first order nature of the phase 
transition for (3a > 11, we also simulated the model at (3a = 1.5, where 
we find it to be strongly first order. We performed 10 5 sweeps at (3 — 1.05, 
and 40,000 sweeps at (3 — 1.02 and 1.0535 (fig. 5-a,b,c). Due to the very 
large tunnelling time (> 10 5 sweeps), we were unable to locate the exact 
critical coupling and verify the first order scaling relation. But the evolution 
curves for both \L\ and P at (3 = 1.02, 1.05 and 1.0535 clearly indicate the 
co-existence of two phases seperated by a large barrier. The intermediate 
runs at (3 = 1.03 (50,000 sweeps) and at (3 = 1.04 (100,000 sweeps) also did 
not show any tunnelling. The data at (3a = 1-1 and 1.5 also indicates that 
the Polyakov loop has a discontinuity exactly across the bulk transition line, 
coincident with the discontinuity in the average value of the plaquette. The 
discontinuity in the values of average Polyakov loop (plaquette energy) at 
Pa = 1-1 and 1.5 is 0.23 (0.04) and 0.48 (0.25) respectively. 

5. DISCUSSION 

From fig. 1 and from the results presented in the previous section, it is 
clear that for Pa — 1-0, the deconfinement phase transition joins the bulk 
transition line, and surprisingly remains apparently coincident for all values 
of Pa > 1-0. As far as, we can tell from our simulations, there is really no 
evidence for two transitions. The discontinuity in both the average plaquette 
and the average Polyakov loop is located at the same P and the latter jumps 
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from an essentially zero value to a rather large value at the same point. 
The fluctuations in both the phases are rather small and are not suggestive 
of any superimposed second order phase transition, although the dominant 
first order nature could make it difficult to see them. The scaling behaviour of 
the finite temperature and bulk transitions is different for iV T , N a — > oo. The 
(3 C for the former should move to infinity, whereas it should remain anchored 
at finite value in the latter case. Therefore the joining together of these two 
transitions is a curious phenomenon, which, if it persists in this limit, would 
lead to a paradox. In the following, we discuss three possible scenarios and 
their implications as one approaches the continuum limit of the theory. 

A] The most conservative possibility is that the joining of the two lines is 
accidental for N T = 4 but they will eventually separate out for large 
enough N T . This scenario is consistent with Bhanot and Creutz's in- 
terpretation and the observed universality of the deconfinement phase 
transition from the equality of its scaling exponents with those of the 
3-d Ising model. To explore this possibility a little further, we simu- 
lated the model at (3a = 11 and 1.5 on a 12 3 x 6 lattice also. Due to 
the large simulation time we could only bracket the critical coupling 
by looking at the peak position of the susceptibility and also the be- 
haviour of the histograms. Our estimated (3 C for /5^=1.1, N T = 6 is 
1.3425 ± 0.0025 which should be compared with the critical coupling 
on 8 3 x 4 lattice, (3 C = 1.3270 ± 0.0008, obtained from our data on the 
susceptibilty peaksQ. Therefore, the shift in (3 C in going from 8 3 x 4 to 
12 3 x 6 lattice at (3a = 11 is definitely less than 0.0188. Moreover the 
transitions still remained coincident and first order. At (3a = 1-5, due 
to the large tunnelling time even on the 8 3 x 4 lattice, we ran a hystere- 
sis cycle instead of runs at individual couplings. The forward hysteresis 
cycle at (3a = 1-5, on 12 3 x 6 shows (3 C < 1.065. On the other hand, 
on the 8 3 x 4 lattice the evolution curves (fig. 5-a,b,c) show that the 
(3 C > 1.02. Therefore the shift in (3 C on going from N T = 4 to iV T = 6 is 
also less than 0.045 here. Recall that at Pa — this shift was observed 



to be 0.13||11||, and if the asymptotic scaling relation were to be valid in 
this region, the expected shift in the critical coupling is 0.15. Of course, 
the scaling relations, eqs. (H-ffl) are known to be strongly violated in 



3 The errorbars given are upper limits based on the couplings of nearby runs which were 
seen to be definitely confined and deconfined. 
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this region. Nevertheless, the shifts are too tiny and at least suggest 
that much larger lattices are required to confirm this scenario if one is 
to move away from the bulk transition and see a separate second order 
deconfinement phase transition. 

B] Taking the coincidence of finite temperature deconfinement phase tran- 
sition and the bulk phase transition more seriously, one has basically 
two extreme alternatives. Either the interpretation of Bhanot and 
Creutz|| of the bulk transition line was incorrect, and one has only 
a deconfinement transition to deal with, or the identification of the de- 
confinement transition at (3a = needs to be investigated more, and 
one has only a bulk transition line on which the order of the transition 
changes from first to second at the point (1.5, 0.9). Both the alterna- 
tives have their own problems and inconsistencies with the published 
results. Thus, the nonzero string tension reported in Ref. || around 
this point suggests a higher order phase transition at that point while 
the results on the growth of susceptibility, reported both in this work, 



and Ref. [10] are indicative of a second order phase transition for 
Pa < 0.9. Indeed, if the entire phase transition line is bulk then all the 
evidence for the universality class of the deconfinement phase transition 
at (3a = will have to be treated as accidental and a priori the global 
Z 2 center symmetry has also no unique role to play in that case. This is 
because for a four dimensional bulk system there is no particular reason 
to choose periodic boundary conditions, from which the Z 2 symmetry 
arises, whereas this is required for the finite temperature interpreta- 
tion. However it should be pointed out that the use of an asymmetric 
lattice will mask the true scaling behavior of a 4-d bulk system, be- 
cause the system will behave as a thin film or layer system which will 
have a three dimensional critical behavior in the limit N a — > oo, N T 



finite [£0] . From this point of view it is not surprising that the system 
exhibits a 3-d critical behavior, typified by the 3-d Ising model, in this 
limit, even if the underlying transition is a 4-d bulk transition which 
would presumably have a completely different critical behavior in the 
limit N T = N a ^ oo. Of course, a crucial prediction of this hypothesis 
is that the entire transition line will persist at finite f3 c as N T — > oo for 
all Pa, and no connection of the confining phase of the pure gluonic 
theory with its asymptotically free phase will be possible. Therefore 
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the continuum limit will be always deconfining, a possibility which has 
been raised previously in several contexts pi]] . 

On the other hand, if the entire line is a deconfinement phase transition 
line, it will move towards the right hand side of the phase diagram as 
N T — > oo. However, our observation of an extremely small shift in (3 C as 
iV T went from 4 to 6 at Pa = 1.1 makes it unlikely to observe this shift 
for any N T used in any numerical simulations up to now. Moreover, the 
large discontinuity in the plaquette and the Polyakov loop then suggests 
the deconfmement phase transition to be of first order at Pa =1.1 and 
1.5. Thus, the universality class seems to change with an apparently 
irrelevant coupling on these lattices. One can then only hope that the 
universality will be restored on larger lattices. 

To summarize, we have studied the finite temperature deconfinement 
phase transition for the extended SU(2) action on N% x iV T lattices with 
N a = 8, 10, 12 and N T = 4 and 6. We found that as Pa increases, the 
deconfinement phase transition moves towards smaller P and appears to join 
the previously known bulk transition line. The universality of the critical 
exponents for the 577(2) deconfinement phase transition thus seems to be 
violated since the finite temperature deconfinement phase transition appears 
to change from being a second order one to a first order transition unless this 
behaviour is purely accidental on small lattices. Future simulations on larger 
lattices could, of course, yield a separation of the bulk transition from the 
deconfinement phase transition, however, no such hints could be obtained in 
the limited variations studied here. It will be interesting to study whether a 
similar phenomenon also takes place for the 577(3) theory. We conclude that 
much larger lattices are required to settle the issues raised in this paper. 
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FIGURE CAPTIONS 



Fig.l The phase diagram of the extended SU(2) lattice gauge theory. The 
solid lines are from simulations done on a 5 4 lattice by Bhanot and Creutz||. 
The broken line is the finite temperature deconfinement phase transition line 
based on the results discussed in the paper. 

Fig.2 Average Polyakov loop {(\L\)) at (a) (3 A =0.5, (b) #4=0.75, (c) (3 A =0.9 
and (d) Pa=1.1 on 8 3 x 4 and 10 3 x 4 lattices. At (3 A = 0.9 and 1.1, results 
on 12 3 x 4 lattice are also shown. The points with error bars are results of 
simulations and the curves are extrapolations by the Ferrenberg-Swendsen 
technique from the run closest to the peak. 

Fig. 3 The susceptibility curves at (a)(3 A =0.5, (b)/?A=0.75, {c)[3 A =0.9 and 
(d)/3 A = 11 on 8 3 x 4, 10 3 x 4 and 12 3 x 4 lattices. The first order (F.O) and 
second order (S.O) predictions, explained in the text, are explicitly shown. 
The higher (lower) predictions are for the 12 3 x 4 ( 10 3 x 4 ) lattice. The 
points and curves are as in Fig. 2. 

Fig.4 The probability density of \L\ at (3 A = 1.1 on 8 3 x 4 (J3 = 1.32635), 
10 3 x 4 ((3 = 1.3265) and 12 3 x A ((3 = 1.32685) lattices. 

Fig. 5 Evolution of \L\ and (P+0.15) on 8 3 x 4 lattice, at (3a = 1.5 and 
(a)/3 = 1.02, (b)(3 = 1.05, (c)/3 = 1.0535. The upper broken lines are cold 
starts and the lower solid lines are hot starts. 



Table 1 



The values of (/3, (3a) at which simulations were performed, f3 c and the 
finite size scaling exponent uj. The expected value for u is 1.97 (3.0) if 
the deconfining phase transion is second order (first order). 
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